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Abstract 

We show (analytically and by numerical simulation) that the zero- 
temperature limit of the distribution of the thermopower S of a one- 
dimensional disordered wire in the localized regime is a Lorentzian, with a 
disorder-independent width of 4-7r 3 A;gT/3eA (where T is the temperature and 
A the mean level spacing). Upon raising the temperature the distribution 
crosses over to an exponential form oc exp (— 2\S\eT / A) . We also consider the 
case of a chaotic quantum dot with two single-channel ballistic point contacts. 
The distribution of S then has a cusp at S = and a tail oc In \S\ for 

large S (with [3 = 1, 2 depending on the presence or absence of time-reversal 
symmetry) . 

I. INTRODUCTION 

Thermo-electric transport properties of conductors probe the energy dependence of the 
scattering processes limiting conduction. At low temperatures and in small (mesoscopic) 
systems, elastic impurity scattering is the dominant scattering process. The energy depen- 
dence of the conductance is then a quantum interference effect .0 The derivative dG/dE of 
the conductance with respect to the Fermi energy is measured by the thermopower S, de- 
fined as the ratio — AV/AT of a (small) voltage and temperature difference applied over the 
sample at zero electrical current. Experimental and theoretical studies of the thermopower 
exist for several mesoscopic devices. One finds a series of sharp peaks in the thermopower 
of quantum point contacts,! aperiodic fluctuations in diffusive conductors,! sawtooth oscil- 
lations in quantum dots in the Coulomb blockade regime,! and Aharonov-Bohm oscillations 
in metal rings.! 

Here we study the statistical distribution of the thermopower in two different systems, 
not considered previously: A disordered wire in the localized regime and a chaotic quantum 
dot with ballistic point contacts. A single transmitted mode is assumed in both cases. In 
the disordered wire, conduction takes place by resonant tunneling through localized states. 
The resonances are very narrow and appear at uncorrelated energies. The distributions 
of the thermopower and the conductance are both broad, but otherwise quite different: 
Instead of the log-normal distribution of the conductance! we find a Lorentzian distribution 
for the thermopower. In the quantum dot, the resonances are correlated and the widths 
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are of the same order as the spacings. The correlations are described by random-matrix 
theory,i'@ under the assumption that the classical dynamics in the dot is chaotic. The 
thermopower distribution in this case follows from the distribution of the time-delay matrix 
found recent ly.@ 

The thermopower (at temperature T and Fermi energy E-p) is given by the Cutler-Mott 
formulai'0 

1 JdE(E-E F )G(E)df/dE 

eT JdEG(E)df/dE ' 1 ' ' 

where G is the zero-temperature conductance and / is the Fermi-Dirac distribution function. 
In the limit T — > Eq. (|1 . 1| ) simplifies to 

„ n^klTdG 
~ 3~eG~dE' 

where G and dG/dE are to be evaluated at E — Ep. We consider mainly the zero- 
temperature limit of the thermopower, by studying the dimensionless quantity 

A dG 

a = ^GdE- (L3) 

Here A is the mean level spacing near the Fermi energy. Since we are dealing with single- 
channel conduction, the conductance is related to the transmission probability T(E) by the 
Landauer formulalM 

G(E) = ^T(E). (1.4) 

The problem of the distribution of the thermopower is therefore a problem of the distribution 
of the logarithmic derivative of the transmission probability. 



II. DISORDERED WIRE 

In this section we study a disordered single-mode wire of length L much greater than the 
mean free path I. This is the localized regime. We compute the thermopower distribution 
in the zero-temperature limit. The analytical theory is tested by comparing with a numer- 
ical simulation. The effect of a finite temperature is considered at the end of the section. 
Electron-electron interactions play an important role in one-dimensional conduction, but we 
do not take these into account here. 



A. Analytical theory 

The localization length £(E) [which is of order I and is defined by lim^oo L^ 1 InT(E) = 
—2/£(E)] and the density of states p(E)_\per unit of length in the limit L — >• oo] are related 
by the Herbert- Jones-Thouless formulate 
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J dE'p(E') In \E - E'\ + constant. 



(2.1) 



The additive constant is energy-independent on the scale of the level spacing. Eq. ( [2 . 1| ) 
follows from the Kramers- Kronig relation between the real and imaginary parts of the wave 
number (the real part determining p, the imaginary part £). Neglecting the width of the 
resonances in the large- L limit, the density of states p(E) = L^ 1 J2i 8(E — Ei) is a sum of 
delta functions, and thus 

LA d 1 A^ 1 . , 

(2.2) 



a 



7i V Ei - E F 



7T dEi(E) 

In the localized regime the energy levels E± are uncorrelated, and we assume that they 
are uniformly distributed in a band of width B around Ep. To obtain the distribution of a, 



P(<t) 



n 



B / 2 dE; 



l-B/2 B 

we first compute the Fourier transform 
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a 




P{k) 



dae ika P(a) 



1 

B 



B/2 
-B/2 



dEe lkA/nE 



1 B/A 



-\k\ 



(2.3) 



(2.4) 



where the limit B/A — > oo is taken in the last step. Inverting the Fourier transform, we 
find that the thermopower distribution is a Lorentzian, 

1/tt 



(J 



1 + a 2 



(2.5) 



The "full width at half maximum" of P(cr) is equal to 2, hence it is equal to 47r 3 /c|T/3eA 
for P(S). This width depends on the length L of the system (through A oc l/L), but it 
does not depend on the mean free path I (as long as I L, so that the system remains in 
the localized regime). 



B. Numerical simulation 

In order to check the analytical theory, we performed a numerical simulation using the 
tight-binding Hamiltonian 

W = -f E OJ+iS- + c h+i) + E Vjcfa. (2.6) 

Z 3 3 

The disordered wire was modeled by a chain of lattice constant a, with a random impurity 
potential Vj at each site drawn from a Gaussian distribution of mean zero and variance u 2 . 
The localization length of the wire is given byE! £ = 2(a/u 2 )(w 2 — i?|). We have chosen 
u = 0.075 w, E F = — 0.55 w, such that £ = 248 a, much smaller than L = 8000 a. From the 
scattering matrix we obtained the conductance via the Landauer formula (|1.4j) , and then the 
(dimensionless) thermopower via Eq. (|1.3j ) (with A = 3.3- 10~ 4 u>). The differentiation with 
respect to energy was done numerically, by repeating the calculation at two closely spaced 
values of Ep. As shown in Fig. [I], the agreement with the analytical result is good without 
any adjustable parameters. 
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FIG. 1. Distribution of the dimensionless thermopower a = (A/2ir)d\nT(E)/dE for a 
one-dimensional wire in the localized regime. The histogram is obtained from a numerical simula- 
tion, for a sample length L = 32.3 £. The dashed curve is the Lorentzian ( |2.5| ), being the analytical 
result for L S> £. The inset shows the algebraic tail of the distribution on a logarithmic scale. The 
thermopower S, in the zero-temperature limit, is related to a by S = — (2^ /3)(k B T / eA)a . 



C. Finite temperatures 



Our derivation of the Lorentzian distribution of the thermopower holds if the temperature 
is so low that k B T is small compared to the typical width 7 of the transmission resonances. 
What if k B T > 7, but still k B T <C A (so that the discreteness of the spectrum remains 
resolved)? We will show that the distribution crosses over to an exponential, but in a highly 
non- uniform way. 

Consider arbitrary 7 and k B T, both <C A. The Cutler-Mott formula (|1.1|) is dominated 
by two contributions, one from a peak in df/dE of width k B T around Ep and one from a 
peak in G(E) of width 70 around Eq. Here 70 and Eq are the width and position of the 
level closest to E F . If \E-p — Eq\ ^> max (k B T, 70), the two peaks do not overlap and one can 
estimate the thermopower as 



eT 



nio(k B T) 2 £ F -£o c -| EF - Eo |/ fcB T 



3(E F - E ) 3 k B T 



x 



7o 



n -1 



+ 



\E F -E \/k B T 



2tt(E f - E ) 2 k B T 



(2.7) 



If k B T <C 70, the first terms in numerator and denominator dominate over the second terms. 
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This is the regime that the Lorentzian distribution ( |2.5|) holds for all S. 

We now turn to the regime k^T > 70 . The first terms dominate if \Ep — E \ 3> 
h^T In fc B T '/ '70 Hence P(S) is a Lorentzian for \S\ <C (^b/c) (In /^bT/^o)" 1 . The logarithm 
of /cbT/7q can be quite large, because the width of the levels is exponentially small in the 
system size, 7 ~ e~ L ^. The Lorentzian persists in an interval larger than its width, provided 
ksT < A (In k&T/j )~ . The second terms in Eq. (|2T?| ) dominate if k^T <C \Ep — E \ <C 
/ceTln /cbT/7 . In this case the thermopower is simply S = (E-p — E )/eT, with exponential 
distribution 

The distribution (|2.8|) follows because the energy levels are uncorrelated, so that the spacing 
\Ef — Eq\ has an exponential distribution with a mean of A/2. 




FIG. 2. Thermopower distribution of a one-dimensional wire in the localized regime at finite 



temperature. The histogram is obtained from Eqs. (1.1) and ( p.lOj ), by numerical integration for 
a set of randomly chosen energy levels Ei, all having the same width 7i = 7 = 10~ 6 A. The 
temperature is feeT/A = 0.01, such that 7 <C k-gT <C A. The distribution follows the Lorentzian 
( |2,5| ) (solid curve) for small and large S, but it follows the exponential fl2.8|) (dashed curve) in an 
intermediate region. 

We conclude that the thermopower distribution for 7 < k^T <C A contains both 
Lorentzian and exponential contributions. The peak region \S\ <C (&b/c) (Ink^T j^)" 1 is 
the Lorentzian (|2.5|). The intermediate region (fca/e) (In /cb^/7) 1 "C \S\ <C (A^B/e) lnksT/^ 
is the exponential (|2~8|). The far tails \S\ 3> ik-Q/e) In k^Tj^ can not be explained by Eq. 
( p.7[ ). With increasing temperature, the Lorentzian peak region shrinks, and ultimately 
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the exponential region starts right at S = 0. This applies to the temperature range 
A (In /c B T/7) _1 < k B T <C A. 

To illustrate these various regimes, we computed P{S) numerically from Eq. ( |1.1| ). We 
took the density of states 

^-'xi (£ _ff; 7?/4 . (2.9) 

so that the conductance according to Eq. ( fOP has the energy dependence 

G{E)^\{[{E-E l f + 1 l/A]' 1 . (2.10) 

The levels Ei were chosen uniformly and independently (mean spacing A), but the fluctu- 
ations of the widths 7; were ignored (7, = 7 for all i). Such fluctuations are irrelevant in 
the low-temperature limit k B T <C 7, but not for 7 < k B T <C A. We believe that ignoring 
fluctuations in 7; should still be a reasonable approximation, because 70 appears only in 
logarithms. The resulting P(S) is plotted in Fig. |2|. We see the expected crossover from 
a Lorentzian to an exponential. The exponential region appears as a plateau. Beyond the 
exponential region, the distribution appears to return to the Lorentzian form. We have no 
explanation for this far tail. 



III. CHAOTIC QUANTUM DOT 

In this section we consider a chaotic quantum dot with single-channel ballistic point 
contacts (see Fig. |3], inset). Because there are no tunnel barriers in the point contacts, the 
effects of the Coulomb blockade are small and here we ignore them altogether. For this 
system, the distribution of dT/dE was computed recently from random-matrix theory.! The 
energy derivative of the transmission probability has the parametrization 

§ = 5fr-W5x^n (3-D 

with independent distributions 

P(c) oc (1 - c 2 )- 1+m , |c| < 1, (3.2) 

P(n,r 2 ) oc \ n - r2 |/3 (riT2) -^ + i) e -(i/n + iM)^/A 5 Ti;T2 > 0> (3 3) 

P(T) oc T~ 1+/3/2 , < T < 1. (3.4) 

The integer (3 equals 1 or 2, depending on whether time-reversal symmetry is present or 
not. The times 7"i,T2 are the eigenvalues of the Wigner-Smith time-delay matrix (see Refs. 
|8] and |1^). Their sum T\ + r 2 is the density of states (multiplied by 2nh). The thermopower 



distribution follows from 



dcP{c) / d n / dr 2 P{T U r 2 ) / dTP{T) 
-1 Jo Jo Jo 



x (n + t 2 )5 a - (A/2nh)c(n - r 2 ) Jl/T - 1 . (3.5) 
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FIG. 3. Distribution of the dimensionless thermopower of a chaotic cavity with two sin- 
gle-channel ballistic point contacts (inset), computed from Eq. (3.5) for the case of broken (/? = 2) 
and unbroken (/? = 1) time-reversal symmetry. 



As in Refs. [8] and [L5| the density of states appears as a weight factor n + t<i in the ensemble 
average fl3.5| ), because the ensemble is generated by uniformly varying the charge on the 
quantum dot rather than its Fermi energy. This is the correct thing to do in the Hartree 
(self-consistent potential) approximation. A more sophisticated treatment of the electron- 
electron interactions (as advocated in Ref. [16]) does not yet exist for this problem. The 
resulting distributions are plotted in Fig. [| The curves have a cusp at a = 0, and asymptotes 
P(a) oc |o-h 1_/3 ln Icrl for \a\ > 1. 



IV. CONCLUSION 



The results we have reported hold for single-channel conductors. The generalization to 
multi-channel conductors is of interest. Multi-channel diffusive conductors were studied in 
Ref. |||. For a chaotic cavity with ballistic point contacts having a large number of modes 
(N modes per point contact), the distribution of the thermopower is Gaussian. The mean 
is zero and the variance is 



VarS 



fc|T 2 7T 6 



(4.1) 



9e 2 iV 4 A 2 /r 

(We have used the results of Ref. [H].) Analogously to universal conductance fluctuations, 
the variance of the thermopower is reduced by a factor of 2 upon breaking time-reversal 
symmetry {(3 — 1 — > /3 — 2). 
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For an iV-mode wire in the localized regime, our derivation of the exponential distribution 
of the thermopower remains valid. This is not true for the Lorentzian distribution. The 
reason is that the Herbert- Jones- Thouless formula for N > 1 relates the density of states 
to the sum of the inverse localization lengths]!! and there is no simple relation between this 
sum and the thermopower. We expect that the tail of the distribution remains quadratic, 
P(S) oc S~ 2 — because of the argument of Sec. [II C[ which is still valid for N > 1. It remains 
a challenge to determine analytically the entire thermopower distribution of a multi-channel 
disordered wire. 
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